Stretching chimeric DNA: a test for the putative S-form 
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Double-stranded DNA 'overstretches' at a pulling force of about 65 pN, increasing in length by a 
factor of 1.7. The nature of the overstretched state is unknown, despite its considerable importance 
for DNA's biological function and technological application. Overstretching is thought by some to 
be a force-induced denaturation, and by others to consist of a transition to an elongated, hybridized 
state called S-DNA. Within a statistical mechanical model we consider the effect upon overstretching 
of extreme sequence heterogeneity. 'Chimeric' sequences possessing halves of markedly different AT 
composition elongate under fixed external conditions via distinct, spatially segregated transitions. 
The corresponding force-extension data display two plateaux at forces whose difference varies with 
pulling rate in a manner that depends qualitatively upon whether the hybridized S-form is accessible. 
This observation implies a test for S-DNA that could be performed in experiment. Our results 
suggest that qualitatively different, spatially segregated conformational transitions can occur at a 
single thermodynamic state within single molecules of DNA. 

PACS numbers: 



INTRODUCTION 

Double-stranded DNA elongates abruptly at a force of 
about 65 pN if it is pulled along its axis [TJ [2] . The result- 
ing 'overstretched' form of the molecule is approximately 
1.7 times longer than helical B-DNA. Overstretching is 
of crucial importance for the biological function of DNA: 
the bacterial protein RecA elongates DNA by a fac- 
tor of 1.5 upon binding [3l [4] [5], a mechanism central 
to homologous recombination and to chromosomal seg- 
regation during cell division [6 j. However, the nature 
of the overstretched state remains a source of consid- 
erable controversy: some think overstretched DNA in 
vitro to be a hybridized form called S-DNA (the 'B- 
to-S' picture) [7J [8] H EE ED E2], while a competing 
picture considers overstretching to signal a conversion 
to unhybrizided single strands (the 'force-melting' pic- 
ture) [T3J OH OH EH [171 [IS]. We therefore lack full 
understanding of the basic mechanical and thermody- 
namic properties of a molecule of central importance to 
biology and of rapidly increasing importance to technol- 
ogy na na no na . 

It is not possible to conclusively validate or rule out 
either picture of overstretching on the basis of existing 
thermodynamic force-extension data. The variation of 
overstretching force with parameters such as tempera- 
ture and salt concentration [T3l [14] US], both of which 
are known to change the melting properties of DNA, 
implies that melting of the double helix plays an im- 
portant role in its elongation. Combining these data 
with the observation that 65 pN of tension provides 
just enough mechanical energy to render B-DNA single- 
stranded or 'molten' [T3] strongly suggests that over- 
stretching is melting, induced by a pulling force. How- 
ever, some authors [8] charge that molten DNA, in the 



sense of two parallel but noninteracting single strands, is 
thermodynamically unstable to a form of DNA in which 
one strand has frayed or 'unpeeled' from the other, and 
no longer bears tension. The latter state may arise when 
melting occurs in the vicinity of 'nicks' in the phosphate 
backbone. Nicks are present in many overstretching ex- 
periments (see, however, Ref. [23]), and are considered to 
relieve the torsional constrant on the structure, so per- 
mitting the two strands to rotate freely when stretched. 
The instability of melting to unpeeling implies that for 
nicked DNA the noninteracting parallel-strand confor- 
mation is not the stable state at 65 pN. However, this 
stable state appears not to be the unpeeled form of the 
molecule, either. Under some conditions overstretched 
DNA possesses a mechanical stiffness far in excess of 
that of single-stranded forms of the molecule [8], imply- 
ing that the stable overstretched form at 65 pN is double- 
stranded, as in the S-DNA picture. This view is bolstered 
by the observation that a transition to single-stranded 
(unpeeled) DNA is sometimes observed only at forces 
well in excess of 100 pN [24[[25]. But these observations 
pose a problem: assuming that overstretching involves 
the inter conversion of two hybridized states, which pre- 
sumably differ principally in their enthalpic rather than 
their entropic properties, how does one rationalize the 
striking temperature dependence of overstretching? 

Just as thermodynamic data fail to distinguish be- 
tween force-melting and B-to-S pictures, so current imag- 
ing techniques cannot resolve the structure of over- 
stretched DNA. Further, atomistic simulations of over- 
stretching [26j [27j [28], while offering valuable insight 
into molecular mechanisms, cannot approach (by orders 
of magnitude) the length and time scales characteristic 
of experiment. We argue, by contrast, that the kinetics 
of force-extension data, combined with predictions from 
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coarse-grained, statistical mechanical models, may offer 
a means of discriminating between the two scenarios. 

In typical experiments, an optical tweezers [29, [30] 01 
atomic force microscope (AFM) is used to overstretch 
a single molecule of DNA at constant rate of exten- 
sion. The loading rate is then reversed, allowing the 
molecule to recover its original length. If this stretching- 
shortening cycle is carried out at low temperature, the 
force-extension traces for each stage superpose, indicat- 
ing a reversible transition. As temperature increases, 
stretching and shortening traces become distinct, signal- 
ing hysteresis. The degree of this hysteresis increases as 
temperature increases [31]. 

Such temperature-dependent kinetics provides a strin- 
gent test of theories of overstretching. Hysteresis is very 
unusual in a system whose extent is macroscopic in only 
one dimension, and whose fluctuations involve local free 
energy differences not much greater than the thermal 
energy, k^T. Well-known causes of hysteresis in phys- 
ical systems include strong local interactions that lead to 
an unfavorable surface tension between coexisting phases 
within, for instance, magnetic materials or liquid-vapor 
phase transitions. Overcoming this surface tension in or- 
der to nucleate and grow domains of the thermodynami- 
cally dominant phase results in sluggish kinetics and hys- 
teresis. However, hysteresis associated with strong inter- 
actions increases in degree as temperature decreases, by 
contrast with overstretching hysteresis. Moreover, sur- 
face tension does not grow with domain size in a quasi- 
one dimensional structure such as DNA. One is led to the 
conclusion that hysteresis in stretching data implies the 
emergence of concerted, long-wavelength correlations. 

Long-wavelength correlations can emerge from the de- 
tachment and re-annealing of strands [TJ [14]. The au- 
thors of Ref. [8 demonstrated that a kinetic model of 
strand separation (unpeeling) displayed hysteresis simi- 
lar in character to that seen in experiment. Building on 
these observations, we introduced a discrete statistical 
mechanical model designed to assess the kinetics associ- 
ated with both the 'B-to-S' and 'force-melting' pictures 
of overstretching [32] . The model is resolved at the level 
of individual basepairs and assumes that DNA may lo- 
cally adopt certain discrete conformational states. It is 
inspired by and borrows features from models of DNA 
undergoing thermal melting [33, 34, 35 and overstretch- 
ing [8] |14j . Our conclusion is that indeed long- wavelength 
correlations of the kind propagated by separating and re- 
annealing strands induce hysteresis. Crucially, in order 
to observe a progression with temperature of the degree 
of hysteresis, we require a progression with temperature 
of the nature of the overstretched state. This progression 
results within our model from a competition between the 
basepairing energy of S-DNA and the entropy liberated 
upon unpeeling. At low temperatures, S-DNA predomi- 
nates as the overstretched form. Our model of the B-to-S 
transition involves only local free energy barriers of mag- 



nitude < 10 &bT, and at pulling rates considered occurs 
in equilibrium. At high temperature, unpeeling occurs, 
whose associated long-wavelength correlations give rise 
to hysteresis. 
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FIG. 1: Overstretching natural and designed sequences, (a) 
Simulated A-DNA (40 kbp, 150 mM NaCl, 21°C, pullrates 
150 and 700 nm/s) displays a B-to-S transition at 65 pN fol- 
lowed at higher forces by a rate-dependent unpeeling. (b,c) 
By contrast, simulations of strongly heterogeneous (chimeric) 
sequences (10 kbp, 150 mM NaCl, 16°C, d = 0, pullrates 
400 and 1500 nm/s) exhibit both transitions in a spatially 
segregated fashion. Corresponding force-extension data dis- 
play double plateaux, whose forces vary strongly with pulling 
rate if plateaux correspond to unpeeling, and are insensitive 
to pulling rate if plateaux signal B-to-S conversion. When S- 
DNA is included in the model, the AT-rich half can unpeel (b, 
lower plateau) and the CG-rich half can convert to S-DNA (b, 
upper plateau); when S-DNA is disallowed both halves unpeel 
(c). Inset: force as a function of temperature at which the 
basepair doublets CG:CG and AT:AT unpeel in preference to 
remaining as B-DNA. These unpeeling transitions compete 
with B-to-S conversion at 65 pN. 

Consequently, we interpret overstretching as a process 
involving competing conformational transitions, with el- 
ements of both pictures required to explain experimen- 
tal results. This view has several qualitiative implica- 
tions for experiment, such as the prediction that hys- 
teresis should be a non-monotonic function of pulling 
rate [32 j. Here we examine the consequences of this 
picture of overstretching in a new context, by explor- 
ing within our model the effect of extreme sequence het- 
erogeneity upon the competition between unhybridized 
and hybridized elongated states. We are motivated by 
the expectation that appropriate sequence design could 
lead, within the same molecule, to distinct transitions 
from B-DNA to two different overstretched states. In 
stretching experiments at room temperature, bacterial 
A-phage DNA (A-DNA) (roughly 50:50 AT:CG content, 
without long-ranged sequence correlations [39]) displays 
a force-extension plateau at about 65 pN, independent of 
pulling rate between 150 and 3000 nm/s. A second tran- 
sition is observed at higher forces, strongly dependent 
upon pulling rate [24] [25] . This behavior is illustrated 
by model calculations in Figure 1, left panel. These ob- 
servations have been rationalized as two, temporally seg- 
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regated transitions, the first an equilibrium B-to-S con- 
version, the second signaling out-of-equilibrium unpeel- 
ing [8j [32] . We demonstrate here that one can dramati- 
cally exaggerate the differences between these transitions 
by designing sequences that permit such changes to occur 
in a spatially segregated fashion. We suggest that testing 
for similar segregation in experiment will reveal if indeed 
there exists an elongated, hybridized state. 



MODEL 

Our model resolves detail at the level of individual 
basepairs. We assume that each basepair, of type AT, 
TA, CG or GC (we assume fully complementary align- 
ment), may instantaneously adopt one of four discrete 
states. These are the helical B-form; two different unhy- 
bridized forms (M, or molten, corresponding to internal 
molten bubbles; and U, or unpeeled, in which only one 
strand is load-bearing); and the putative S-form. The 
model takes as its input coarse-grained free energetic 
properties of each state, and via a Monte Carlo algo- 
rithm offers kinetic predictions on length- and timescales 
characteristic of experiments. Full details are given in 
Ref. [32], with two modifications considered in this paper. 
The first accounts for a nearest-neighbor dependence of 
basepairing-stacking energies [35 , while the second ac- 
counts for the sequence-dependence of the heat capacity 
of melting p~4[ [15]. We assume that specific heats of 
melting of AT and CG basepairs are zero at (i.e. Tay- 
lor expanded about), respectively, the (salt-dependent) 
melting temperatures of AT- and CG-DNA. This gives a 
model in which the stabilities of AT and CG basepairs 
are more similar at lower temperatures or higher salt con- 
ditions than under the converse conditions. While such 
details matter greatly when considering changes in over- 
stretching behavior with temperature, they play essen- 
tially no role in determining changes with pulling rate. 
We evolve our model according to a dynamic protocol 
designed to mimic optical trap or AFM pulling experi- 
ments. We increment at constant rate the position of one 
end of the molecule, allowing the position of the molecule 
end tethered to an imaginary optical trap or cantilever 
to fluctuate. We calculate the resulting tension subject 
to the constraint of mechanical equilibrium. Basepair 
fluctuations are assumed to occur on a timescale derived 
from fluorescence correlation spectroscopy [36j [37l [38] : 
we discuss this choice in Appendix B. 



RESULTS 

We consider simulations of 'chimeric' DNA fragments 
whose left halves consist of basepairs having chemi- 
cal composition AT or TA with respective probabilities 
(1 — d)/2, and CG or GC with respective probabilities 
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FIG. 2: Simulated force-extension data for strongly heteroge- 
neous sequences bear a kinetic 'fingerprint' of the nature of the 
spatially segregated overstretching mechanisms. The pulling 
rate-dependence of A/ P i a t for 10 kbp chimeric sequences at 
16° C and 150 mM NaCl (d — [each data point is the mean 
of 5 simulations]; and d — 0.02 [each data point is the mean 
of 2 simulations]) is qualitatively altered by the presence of 
an elongated, hybridized state. Plateaux height difference if 
S-DNA is disallowed (a) is a weak function of pulling rate, but 
is a strong function of pulling rate if S-DNA is included in the 
model (b). Absolute plateaux heights border the shaded ar- 
eas (main figure, right axis). Inset: A/ P i a t at 1000 nm/s as a 
function of sequence disorder d within the B-to-S model. The 
double plateaux structure (corresponding to nonzero A/ P i at ) 
is visible for d < 0.15. 



d/2. We call d the 'sequence disorder' parameter. The 
right halves have a corresponding structure with the re- 
placement d — > 1 — d. Thus d = indicates a perfectly 
segregated sequence, while d = 1/2 corresponds to a com- 
pletely random sequence. We place one nick at the left 
extremity of the AT-rich half. The nick permits unpeel- 
ing, which at the pulling rates and temperatures we con- 
sider occurs in preference to the internal molten bubble 
configuration (M) . We perform calculations on sequences 
of lengths from 10 kilobasepairs to 200 basepairs. 

In Figure 1 (b,c) we show that in simulations modeling 
stretching near room temperature these constructs dis- 
play two distinct force-extension plateaux, corresponding 
to unpeeling of the AT-rich half (b,c, lower plateaux) and 
one of two fates for the CG-rich half. When we permit 
the model to access the S-state, the CG half can elon- 
gate by way of the B-to-S transition (b, upper plateau). 
When S-DNA is not included in the model, the CG half 
unpeels, albeit at a higher force than does the AT-rich 
half (c, upper plateau). Plateaux resulting from un- 
peeling, which occurs out of equilibrium, display strong 
pulling rate dependence [8j [24j [32]; plateaux signalling 
equilibrium B-to-S conversion do not. Consequently, the 
difference in 'height' (force) between double plateaux, 
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A /plat ? varies with pulling rate in a manner that depends 
qualitatively upon whether the molecule may access the 
hybridized S-state. When S-DNA is accessible, A/ P i at 
varies strongly with pulling rate. When S-DNA is in- 
accessible we find that the difference in height between 
unpeeling plateaux is almost insensitive to pulling rate, 
because the mechanism of strand separation is similarly 
hysteretic for both halves. We note that at higher forces, 
within the B-to-S model, one can observe a three-stage 
transition within chimeric sequences (Appendix A). We 
note also that at temperatures high enough that S-DNA 
is unstable at about 65 pN to unpeeling of CG basepairs, 
the B-to-S model describes a variation of A/ P i at with 
pulling rate similar to that of the force-melting model 
(data not shown). 




FIG. 3: Microscopic configurations as a function of time 
from simulations of 300-basepair fragments of A-DNA (a) and 
chimeric DNA (b,c) at 21°C and 150 mM NaCl, for pulling 
rate 1000 nm/s. White indicates B-DNA, red indicates un- 
hybridized DNA, and blue indicates S-DNA. The spatially 
distinct transitions observed within chimeric sequences are of 
different character depending on whether S-DNA is assumed 
to exist (b) or not (c). The spatial scale indicates molecu- 
lar length, with the extension modeling the departure of the 
'trapped' bead from the trap center reduced by a factor of 10 
for clarity. The bottom (top) halves of chimeric sequences are 
AT (CG) rich. The maximum force attained is 90 pN. 

In Figure [2] we quantify these observations by plotting 
A/ P iat versus pulling rate with S-DNA included in the 
model (a) or not (b). The behavior of A/ p i at is qualita- 
tively and strikingly different within the two scenarios, 



and constitutes the key result of this paper. The cor- 
responding pulling experiments would reveal, firstly, if 
force-extension data for chimeric sequences exhibit dou- 
ble plateaux, and secondly, if so, how these plateaux vary 
with pulling rate. We argue that pronounced variation 
with pulling rate of A/ P i at would imply overstretching 
via two distinct mechanisms, suggesting that elongation 
is possible via both unhybridized and hybridized states. 

In Figure [3] we show the time-dependent microscopic 
configurations that underlie these force-extension data, 
comparing A-DNA (a) with chimeric DNA (b,c). The 
spatially distinct transitions observed within strongly 
heterogeneous sequences are of different character de- 
pending on whether S-DNA is assumed to exist (b) or 
not (c) . The microscopic signature of unpeeling is a long- 
wavelength drift of the domain wall separating hybridized 
and unhybridized conformations. 

Our results indicate that clean segregation between 
overstretching transitions (signaled by distinct double 
plateaux in force-extension data) occurs for molecules 
as short as 200 bp (data not shown), and for sequence 
disorder parameters d as large as d « 0.15. For the B- 
to-S model, these segregated transitions are different in 
nature (unpeeling versus B-to-S conversion) for temper- 
atures T < 25°C at 150 mM NaCl (and for higher tem- 
peratures at higher salt concentrations). These results 
suggest that DNA-elongating proteins, such as RecA, 
could induce within DNA stretching transitions whose 
character differs with sequence composition, conferring 
upon the elongated Rec A-DNA complex [3jH] a sequence- 
dependent elasticity. Our results also suggest that pre- 
cise melting of localized regions of DNA can be effected 
by stretching molecules whose sequences are appropri- 
ately designed. Such precision de- hybridization would 
facilitate basepair doping with ligands, allowing one to 
modify the conductance properties of the molecule in a 
spatially heterogeneous fashion. Control of conductance 
would enhance DNA's usefulness as a nanowire in molec- 
ular electronics devices. 

We expect that in experiment the effects of secondary 
structure will play only a minor role in the scenarios we 
have discussed. While it is likely that hairpins, for in- 
stance, will form in separated strands that are very AT- 
or CG-rich, and while indeed hairpins give rise to addi- 
tional plateaux at - 9 pN (AT) and ~ 20 pN (CG) upon 
repeated extension [24] , we argue that hairpin formation 
will not influence strongly the heights of the plateaux 
measured upon first extension of the DNA. Plateaux 
heights will be influenced if the dynamics of the unpeeled 
'front' is restricted by the formation of a hairpin in the 
non tension-bearing strand. This would be a concern 
were the front to move diffusively, in which case its re- 
verse movement might be blocked by a hairpin formed in 
its wake. However, at the pulling rates we consider the 
motion of the front is super-diffusive, and its dynamics is 
not strongly affected by potential blockages in its wake. 
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FIG. 4: Temperature dependence of plateaux height differ- 
ence A /piat from simulations of 10 kbp chimeric sequences at 
150 mM NaCl, for d = [each data point represents mean 
of 2 simulations] is not a discriminating test of the two pic- 
tures of overstretching, (a) When S-DNA is assumed not 
to exist the temperature variation of A/ P i a t is governed by 
the temperature dependence of the specific heat of melting of 
AT- and CG basepairs. This variation is weak when specific 
heats are taken to be zero at (i.e. Taylor expanded about) 
the melting temperature of A-DNA (dashed line), and strong 
when specific heats are taken to be zero at respective AT- 
and CG melting temperatures, a more plausible model (solid 
line), (b) When S-DNA is included in the model the AT-rich 
half unpeels at all temperatures considered, while the CG-rich 
half unpeels at high temperature and elongates by way of the 
B-to-S transition at low temperature. Despite this change 
in mechanism upon allowing S-DNA, the resulting plateaux 
height difference varies with temperature in a manner that is 
qualitatively similar to the solid line in the top panel. 



ing test for S-DNA is kinetic in nature. 

We have shown that within a statistical mechanical 
model of DNA overstretching one can induce spatial 
segregation of competing transitions through appropri- 
ate sequence design. The force-extension signatures of 
these distinct transitions vary with pulling rate in a man- 
ner that depends qualitatively on the accessibility of an 
elongated, hybridized state. We propose that the corre- 
sponding experiments would provide a means of assessing 
whether indeed such a state exists. Our results also sug- 
gest that de-hybridization of specific locations within a 
DNA molecule can be effected by subjecting appropri- 
ately designed sequences to external force. Such preci- 
sion melting would permit, for instance, doping of DNA 
at specific locations, allowing fine control of its conduc- 
tance properties. 
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APPENDIX A: MULTI-STAGE ELONGATION OF 
CHIMERIC SEQUENCES 



We end with the observation that varying temperature 
provides a much less clear assessment of the nature of dis- 
tinct overstretching transitions than does varying pulling 
rate. In our simulations, force-extension data for strongly 
heterogeneous sequences display double plateaux, but the 
variation of the heights of these plateaux with tempera- 
ture does not indicate clearly the nature of the respon- 
sible overstretching mechanism. We demonstrate this 
point in Figure |1J Accounting for the sequence depen- 
dence of specific heats of melting, the temperature de- 
pendence of the plateaux height difference is not changed 
qualitatively by allowing or suppressing the S-state. To 
make predictions on the basis of the slight difference that 
does exist would require precise knowledge of the se- 
quence dependence of the specific heat of melting [T5] . 
Indeed, the variation with T of the height of the unpeel- 
ing plateaux (accessible even within the B-to-S model 
at sufficiently high temperature) should allow one to use 
stretching experiments to measure the specific heat of 
melting per basepair. We conclude that the most strik- 



Simulations performed on chimeric sequences in Figure 
1 of the main text (panels b,c) used a maximum pulling 
force of 120 pN, and reveal two distinct plateaux. At 
higher forces we observe within the B-to-S model a third 
plateaux corresponding to a pulling rate-dependent un- 
peeling of the CG-rich half. A similar unpeeling is seen 
when stretching A-DNA to high forces (see Fig 1(a) and 
Refs. [HI [241 [2S]) We illustrate the multi-stage elon- 
gation of chimeric sequences in Figure [5] (a): the low- 
est plateau represents unpeeling of the AT-rich half; the 
middle plateau signals elongation of the CG-rich half by 
way of the B-to-S transition; and the highest plateau 
indicates the subsequent unpeeling of the CG-rich half. 
At pulling rates considered, transitions signaling unpeel- 
ing are pulling rate-dependent, while B-to-S conversion 
is not. In Figure 1(b) we show for comparison calcula- 
tions performed with S-DNA not included in the model. 
In this case the two plateaux signal unpeeling of AT- and 
CG-rich halves; there is no further transition at higher 
force. 

In Figure [6] we show the microscopic dynamics that 
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FIG. 5: Overstretching designed sequences to high forces. 
Simulated chimeric DNA (10 kbp, 150 mM NaCl, 16°C, d = 
0, pullrates 400, 1000 and 1500 nm/s) shows a three-stage 
transition if S-DNA is included in the model (a). The lowest 
plateau signals unpeeling of the AT rich half, while subsequent 
plateaux indicate elongation of the CG-rich half by way of the 
B-to-S transition, followed at higher force by unpeeling of the 
same half. If S-DNA is not included in the model (b) then 
only two transitions are seen. These correspond to unpeeling 
of AT- and CG-rich halves, respectively. 



chimeric, model with S-DNA 




FIG. 6: Microscopic configurations as a function of time 
from simulations of 300-basepair fragments of chimeric DNA 
at 21°C and 150 mM NaCl, for pulling rate 1000 nm/s. S- 
DNA is included in the model. White indicates B-DNA, red 
indicates unhybridized DNA, and blue indicates S-DNA. The 
three-stage elongation evident in the force-extension data of 
Figure [5] (a) is seen clearly as a unpeeling of the AT-rich 
(bottom) half of the molecule, followed by a conversion of the 
CG-rich (upper) half first to S-DNA and then to unpeeled 
DNA. The spatial scale indicates molecular length, with the 
extension modeling the departure of the 'trapped' bead from 
the trap center reduced by a factor of 10 for clarity. The 
maximum force attained is 200 pN. 



underlie the three-stage elongation shown in Figure^] (a). 



fluctuations identified by NMR experiments. Here we il- 
lustrate the effect of changing r by determining analyti- 
cally, as a function of model parameters, the approximate 
force at which complete unpeeling occurs when stretch- 
ing once-nicked DNA. We will compare this estimate with 
the A-DNA overstretching data of Refs. [24] [25]. These 
data display a pulling rate-independent '65 pN' plateau 
followed at higher force by a rate-dependent transition 
interpreted as an unpeeling. 

We consider within our model the drift of the unpeeled 
'front' separating hybridized and unhybridized regions, 
for a sequence of N basepairs bearing a nick at one 
end. For simplicity we will assume sequence homogene- 
ity, and model the force-extension profile of DNA using 
the piecewise linear fit of Cocco et al. [8 j. This fit as- 
sumes a linear overstretching plateau between 62 and 68 
pN of gradient m x = (68 — 62) /(0. 58 — 0.34) pN/(nm- TV), 
and a gradient for forces above 68 pN of magnitude 
?7i2 = 1600/0.34 pN/(nm • N). We shall assume a con- 
stant rate of molecular extension, L = v$t. The mas- 
ter equation derived from our model for the position 
with time t of the unpeeled front, n{t) (0 < n < A/"), 
may be manipulated to yield the front's drift veloc- 
ity, h(t) = T (W+(f) - W-(f)). Here the W ± (f) = 
[1 + exp(±(3[Ag - Sw(f)})}- 1 are Glauber rates. The 
term Ag is the (sequence-averaged) free energy differ- 
ence between hybrizided and unhybridized DNA at zero 
force, and Sw(f) models the free energy of extension of 
the unhybridized phase relative to the hybrdized phase. 
We shall assume that unpeeling occurs only over S-form 
DNA, i.e. that the B-to-S conversion just pre-empts un- 
peeling. For forces large enough that unpeeling is favor- 
able we have 
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APPENDIX B: CHOICE OF MODEL TIMESCALE 

The pulling rate-dependence of unpeeling within our 
model is governed by the imposed pulling timescale, 
controlled by the pulling speed Vq, and the fundamen- 
tal timescale Tq 1 on which basepairs change state. In 
Ref. [32 we argue, by comparison with experimental 
data [24] [25] , that a fluorescence- derived time in the mi- 
crosecond range (r^ 1 =28 /is [37]) is a more appropriate 
choice for this fundamental timescale than the 10-100 ns 



FIG. 7: Approximate unpeeling force f u (derived from Equa- 
tion (|3|) versus pulling speed vq for sequence- averaged A- 
DNA at a range of model timescales Tq 1 . We obtain rea- 
sonable comparison with experiment [24] for timescales in 
the microsecond range, characteristic of fluorescence measure- 
ments [37] . 

dn(t) _ r sinh /3 fAL 
dt ~ °l + cosh/3/Ai' 

where AL « 0.02 nm is the difference in length per base- 
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pair between S- and U-DNA, which we assume to be 
independent of force (a reasonable approximation near 
and just above the overstretching force). The variable 
f = f — / eq- is the force in excess of the S-U equi- 
librium force / eq ., which within our model is approxi- 
mately 62 pN at 150 mM NaCl and 21° C. We assume 
that n(f = 0) = 0. We change variables in Equation |l]) 
from t to L = t/vo, and integrate along the piecewise- 
linear force-extension profile to obtain an expression for 
n in terms of /: 
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^) d? -i In 1 1 + cosh fiALf' I . 
df'J df> 
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This expression may be manipulated to yield 

K[C{h)-C{Q))+C{h) (/>A) 



C(f) 



( k 2 1 n , v ^ v j | ; i ^v./iy v./ — ./ i > /o\ 

1 fcj^n + CtO) (0 < / < A) 



We have defined C(f) = In |l+cosh/?AL/|; h = (68-62) 
pN = 6 pN; ki(v ) = r /(J3ALvomi), for i = 1,2; 
and i\"(v ) = fci/Afe- Recall that /? « 1/(4.1 pNnm). 
Setting n = iV in Equation ([3| gives the excess force 
fu = fu — /eq. a t which complete unpeeling occurs. For 
low pulling rates this excess force scales as /3ALf u w 

1/2 

2 (/^AL^^iAT" 1 ) 7 , while for very large pulling rates 
we have /3ALf u ~ /3ALvo1712NTq 1 . We plot the unpeel- 
ing force f u = f u + / eq . against pulling rate in Figure M 
for different choices of the fundamental timescale Tq . 
We include in this plot data extracted from Figure 2 of 
Ref. [24]. We expect that our estimate of f u (vo) should 
be accurate to within an order of magnitude; we obtain 
favorable comparison with experiment for Tq 1 in the mi- 
crosecond range. The equation of the line in Figure [7] 
calculated with Tq 1 = 5 /is is 



/ = 62 



205 cosh" 1 (1.915 e°- 1148v - l) (/ > 68 pN) 
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0. 0000563i; 



(/<68pN) 
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